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Abstract: 

The  Sanchez-Lacombe  equation  of  state  was  used  to  model  the  sorption  of  high- 
pressure  gases  into  solid,  amorphous  polymers  or  molten  polymers.  One  adjustable 
parameter  per  binary  pair,  812,  was  used  in  the  mixing  brie  to  correct  the  deviation  of  the 
characteristic  pressure  of  the  mixture,  P12*»  from  the  geometric  mean.  The  values  of  512 
which  gave  the  best  fit  of  the  available  literature  data  for  the  carbon  dioxide-polymethyl 
methacrylate,  carbon  dioxide-silicone  rubber,  ethylene-low  density  polyethylene,  methane- 
polyisobutylene,  methane-low  density  and  high  density  polyethylene,  and  methane- 
polystyrene  systems  ranged  from  0.020  to  0.194,  except  for  the  methane-low  density 
system  which  had  a  negative  value  of  -0.183.  In  all  cases,  the  calculated  sorption 
isotherms  were  in  reasonably  good  agreement  with  the  experimental  data.  The  resultant 
swelling  of  polymethyl  methacrylate  and  silicone  rubber  in  the  presence  of  high-pressure 
carbon  dioxide  was  also  well  represented  by  the  Sanchez-Lacombe  equation  of  state. 
Because  the  Sanchez-Lacombe  equation  of  state  is  based  on  lattice-fluid  theory,  the 
sorption  calculations  are  limited  to  polymers  which  are  non-crystalline,  not  cross-linked  or 
slightly  cross-linked,  above  their  glass  transition  temperature,  or  above  their  melting 
temperature.  All  of  the  amorphous  polymers  considered  in  this  study  were  above  their 
glass  transition  temperature  or  their  effective  glass  transition  temperature  as  predicted  from 
a  theoretical  relation  presented  by  Chow. 
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Introduction: 

The  solubility  of  a  gas  in  an  amorphous  or  molten  polymer  is  an  important 
consideration  in  membrane  and  polymer  processes.  For  instance,  the  efficacy  of  a 
membrane  used  for  separating  a  gas  mixture  is  dependent  on  the  solubility  of  the  various 
species  in  the  membrane  (Koros,  1985).  Although  high-pressure  gradients  across  the 
membrane  can  increase  the  solubility  of  one  or  more  of  the  components  from  the  gas 
mixture,  a  very  high  pressure  gradient  can  cause  the  membrane  to  swell  and  the  selectivity 
to  decrease  (Koros,  1985).  An  estimate  of  the  solubility  of  gases  in  the  membrane  over 
large  ranges  of  pressures  could  be  used  to  define  the  operating  boundaries  over  which  the 
membrane  would  effectively  separate  the  constituents  in  a  gas  mixture. 

Swelling  a  solid  polymer  matrix  with  a  high-pressure  gas  can,  however,  aid  in  the 
deposition  of  temperature-sensitive  materials  into  the  polymer.  The  supercritical  fluid 
swells  the  polymer  (Liau  and  McHugh,  1985)  and,  thus,  allows  the  substance  to  migrate 
into  the  polymer  matrix.  Sand  (1986)  has  shown  that  substances  such  as  fragrances,  pest 
control  agents,  and  pharmaceutical  drugs  can  be  impregnated  in  a  solid  polymer  if  the 
polymer  is  exposed  to  a  supercritical  fluid  during  the  impregnation  process.  Swelling  the 
polymer  effectively  increases  the  diffusion  coefficient  of  the  heavy  dopant  by  several  orders 
of  magnitude  and,  thus,  allows  it  to  be  transported  into  the  polymer  within  a  reasonable 
time.  When  the  system  is  depressurized,  the  supercritical  fluid  is  removed  from  the 
polymer  and  the  heavy  dopant,  which  is  trapped  in  the  matrix,  slowly  diffuses  out  of  the 
solid  at  a  rate  which  is  now  orders  of  magnitude  slower  than  the  rate  at  which  it  was  put 
into  the  polymer.  This  technique  could  lead  to  the  formation  of  novel  controlled-release 
devices  since  very  high  molecular  weight  constituents  could  be  impregnated  into  a  swollen 
polymer  matrix  at  operating  temperatures  low  enough  to  avoid  thermal  degradation  of  the 
substituents  if  a  supercritical  fluid  solvent  is  used  as  the  swelling  agent. 
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Of  course  polymers  can  be  effectively  stripped  of  low  molecular  weight  impurities 
including  entrapped  solvent,  residual  catalyst,  or  low  molecular  weight  oligomers  by 
contacting  with  a  high-pressure  gas  or  a  supercritical  fluid.  In  this  instance,  the 
supercritical  fluid  would  swell  the  polymer,  diffuse  into  the  polymer  matrix  and  dissolve 
the  impurity,  and  then  diffuse  out  of  the  swollen  matrix  thus  removing  the  impurity.  This 
supercritical  fluid  extraction  process  is  somewhat  analogous  to  those  used  for  the 
decaffeination  of  coffee  beans  (McHugh  and  Krukonis,  1986)  and  the  removal  of  low 
molecular  weight  constituents  from  solid  coal  particles  (Whitehead  and  Williams,  1975). 

For  any  polymer-gas  process  which  operates  at  high  pressures,  it  is  often  necessary 
to  be  able  to  predict  the  solubility  of  the  gas  in  the  polymer  and  the  subsequent  swelling  of 
the  polymer.  Fleming  and  Koros  (1986)  recommend  using  the  Flory-Huggins  equation  to 
calculate  the  sorption  of  high-pressure  gases  in  rubbery  polymers,  such  as  carbon  dioxide 
in  silicone  rubber.  At  pressures  as  high  as  1000  psia,  they  calculate  the  activity  of  carbon 
dioxide  as  the  ratio  of  system  pressure  to  the  vapor  pressure  of  the  gas.  They  recognize 
that  at  these  high  pressures  the  activity  should  be  the  ratio  of  the  fugacity  of  the  gas  to  the 
standard  state  fugacity.  When  the  ratio  of  system  pressure  to  the  vapor  pressure  of  the  gas 
is  used,  it  is  necessary  to  determine  a  hypothetical  vapor  pressure  of  the  gas  if  the  gas  is  at 
conditions  above  its  critical  properties.  The  focus  of  this  paper  is  to  report  on  the  use  of  the 
Sanchez-Lacombe  equation  of  state  (1976)  for  modeling  the  solubility  of  high-pressure 
gases  in  solid,  amorphous  polymers  or  molten  polymers  and  the  subsequent  swelling  of  the 
polymer.  By  using  an  equation  of  state,  we  avoid  the  need  to  calculate  a  hypothetical  vapor 
pressure  and  the  activity  of  the  high-pressure  gas  can  be  calculated  in  a  thermodynamically 
rigorous  manner.  The  calculational  methods  that  are  described  are  applicable  to  polymers 
which  are  non-crystalline,  not  cross-linked  or  slightly  cross-linked,  above  their  respective 
glass  transition  temperatures  (Tg),  or  above  their  melting  temperature.  The  systems 
modeled  in  this  study  include  carbon  dioxide-polymethyl  methacrylate  (PMMA)  (Liau  and 
McHugh,  1985),  carbon  dioxide-silicone  rubber  (Fleming  and  Koros,  1986),  ethylene- 


4 


low  density  polyethylene  (LDPE)  (Cheng  and  Bonner,  1977),  methane-polyisobutylene 
(PIB)  (Lundberg,  et  al.,  1969),  methane-LDPE  (Lundberg,  et  al.,  1969),  methane-high 
density  polyethylene  (HDPE)  (Lundberg,  et  al.,  1969),  and  methane-polystyrene  (PS) 
(Lundberg,  et  al.,  1969). 

Modeling: 

Bonner  (1977)  presents  a  review  of  the  available  literature  information  on  the 
solubility  of  supercritical  fluids  in  polymers.  He  also  presents  the  predictive  techniques 
available  for  calculating  gas  solubilities  in  polymers.  Our  methods  closely  follows  those 
outlined  by  Bonner  (1977). 

* 

Our  approach  is  to  model  the  polymer-gas  system  using  the  following  assumptions: 

1.  the  solubility  data  of  the  gas  in  the  polymer  represents  equilibrium  data; 

2.  the  solubility  of  the  polymer  in  the  high-pressure  gas  phase  is  essentially  zero; 
and 

3.  the  amoiphous,  solid  polymer  above  its  Tg  can  be  modeled  as  a  liquid. 

The  first  assumption  will  only  be  true  if  the  polymer  is  above  its  glass  transition 
temperature  and  if  sufficient  time  has  been  allowed  for  the  system  to  reach  equilibrium.  A 
number  of  studies  have  shown  that  if  an  amoiphous  polymer  sorbs  gas  to  any  great  amount 
its  glass  transition  temperature  can  be  substantially  lowered  (Chow,  1980;  Wang,  Kramer, 
and  Sachse,  1982;  Liau  and  McHugh,  1985;  Koros,  1985;  Wissinger  and  Paulaitis,  1986). 
For  the  polymers  considered  in  this  paper,  either  the  system  temperature  is  already  above 
the  polymer's  Tg,  or  the  system  temperature  is  above  the  polymer's  effective  Tg  as  gas  is 
solubilized  in  the  polymer. 

The  second  assumption,  that  the  polymer  is  essentially  insoluble  in  the  high- 
pressure  gas,  is  satisfied  if  the  polymer  has  a  very  high  molecular  weight  and  if  the 
molecular  weight  distribution  is  reasonably  narrow.  The  small  amount  of  data  available  in 
the  literature  on  the  solubility  of  polymers  in  supercritical  fluids  suggests  that  if  the 
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molecular  weight  of  the  polymer  is  above  approximately  8,000,  then  the  solubility  will  be 
very  low  (Ehrlich  and  Graham,  1960;  Heller,  et  al.,  1983;  McHugh  and  Krukonis,  1986; 
Dhalewadikar,  McHugi ,  and  Guckes,  1987);  hence,  we  will  assume  this  assumption  to  be 
valid  in  our  analysis. 

The  equation  of  state  of  Sanchez-Lacombe  (1976)  is  based  on  the  lattice-fluid 
theory  which  assumes  that  the  polymer  has  a  flexible  liquid  structure.  This  equation  of 
state  is 

p2+P+f{ln{l-‘p}  +{l-i}p}  =0  (1) 

where  P,T,  and  p  are  the  reduced  pressure,  temperature,  and  density,  respectively,  and  r 
represents  the  number  of  lattice  sites  occupied  by  a  molecule.  The  characteristic  reduction 
parameters  used  for  these  reduced  properties  are  obtained  from  pure-component,  P-V-T 
data  as  described  in  the  following  paragraphs. 

At  equilibrium,  the  chemical  potential  of  the  gas  in  the  pure  gas  phase  is  equal  to  the 
chemical  potential  of  the  gas  dissolved  in  the  polymer 

H°<T,P)  =  ^(T,P,<t>)  (2) 

where  the  superscripts  G  and  P  represent  the  gas  and  polymer  phases,  respectively, 
subscript  1  represents  the  light  gas,  and  represents  the  volume  fraction  of  component  i. 

Using  the  Sanchez-Lacombe  equation  of  state,  the  chemical  potential  of  the  gas  dissolved  in 
the  polymer  can  be  obtained 

it  =  KT{ln  v(l-  rK+V  Xi#  +r*Kr  {.jl  .iil  ♦ 

i  i  r2zi1'4 

?{{l-P}ln{l-p}  ♦ilnpr}} 


(3) 
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where  ij  is  the  number  of  lattice  sites  occupied  by  molecule  1  and  Xj  is  a  mixture  parameter 
which  is  expressed  as 


X 


l 


AP*v1* 

kT 


(5) 


where  AP*  is  subsequently  defined.  The  close-packed,  volume  fraction  of  component  i, 
<j>i,  is  given  by 

TOj/p.* 

«>.  = - 2 -  (4) 

1  mj/p^+mj  Ip* 

where  ra*  is  the  weight  fraction  and  pj*  is  the  characteristic  density.  Although  Sanchez  and 
Lacombe  (1976)  distinguish  between  r^,  the  number  of  lattice  sites  occupied  by  a  molecule 
in  a  pure  state,  and  r2,  the  number  of  lattice  sites  occupied  by  a  molecule  in  a  mixture,  these 

quantities  are  assumed  to  be  equal  in  our  analysis. 

The  chemical  potential  of  the  gas  in  the  pure  gas  phase  is  obtained  by  setting  <j>2 

equal  to  zero  in  Equation  (3).  The  amount  of  gas  sorbed  in  the  polymer  can  be  determined 
by  solving  equation  (2)  with  pure  component  characteristic  parameters  and  with  the 
appropriate  mixing  rules  for  the  mixture  characteristic  parameters. 

The  pure  component  parameters  used  in  this  study  are  listed  in  Table  1.  The 
characteristic  parameters,  p*,  T*,  and  P*,  for  PMMA,  PDMS,  LDPE,  HDPE,  PIB,  PS, 
and  methane  are  obtained  directly  from  Sanchez  and  Lacombe  (1978).  The  characteristic 
parameters  for  PDMS  are  used  as  the  characteristic  parameters  for  the  silicone  rubber,  as 
explained  later. 

As  shown  in  Table  1,  the  characteristic  parameters  for  carbon  dioxide  determined 
from  a  fit  of  the  vapor-liquid  equilibrium  curve,  did  not  follow  the  logical  succession  of 
pure-component  parameters  of  other  light  gases  given  by  Sanchez  and  Lacombe  (1978). 
Thus,  the  parameters  for  carbon  dioxide  used  in  the  calculations  presented  in  this  paper 
were  obtained  by  fitting  data  on  the  liquid  molar  volume,  the  isothermal  compressibility, 
and  the  thermal  expansion  coefficient  of  carbon  dioxide  at  50,  60,  and  75  °C  and  150  atm 
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(Michels  and  Michels,  1937)  to  the  equation  of  state.  These  pressure-temperature  (P-T) 
conditions  were  chosen  so  the  parameters  would  be  based  on  the  P-T  space  in  which  most 
of  the  experimental  polymer-carbon  dioxide  data  were  obtained.  Although  the  parameters 
for  carbon  dioxide  found  in  this  manner  are  more  consistent  with  those  recommended  by 
Sanchez  and  Lacombe  for  similar  light  gases,  we  only  obtain  a  fair  representation  of  the 
vapor-liquid  equilibrium  curve  and  the  critical  point  for  carbon  dioxide  (Tccalc=42.9  °C, 
Pcca,c-96.5  atm;  Tc=31.05  °C,  Pc=72.8  atm). 

Contrary  to  the  carbon  dioxide  case,  the  parameters  for  ethylene  as  determined  from 
a  fit  of  the  vapor-liquid  equilibrium  line  did  follow  the  trends  of  the  parameters  tabulated  by 
Sanchez  and  Lacombe  (1976).  However,  when  these  parameters  are  used  to  calculate 
sorption  isotherms,  it  is  necessary  to  use  a  negative  value  for  the  mixture  parameter,  as 
described  in  the  following  paragraphs.  Therefore,  the  pure  component  parameters  for 
ethylene  are  determined  by  fitting  the  liquid  molar  volume,  the  isothermal  compressibility, 
and  the  thermal  expansion  coefficient  at  135  and  140  °C  and  50  atm  (Angus,  Armstrong, 
and  de  Reuck,  1972)  to  the  equation  of  state.  These  P-T  conditions  are  close  to  those  for 
most  of  the  ethylene-LDPE  experimental  data.  The  parameters  for  ethylene  found  in  this 
manner  are  consistent  with  those  recommended  by  Sanchez  and  Lacombe  for  similar  light 
gases,  although  they  result  in  only  a  fair  representation  of  the  vapor-liquid  equilibrium 
curve  and  the  critical  point  for  ethylene  (Tccalc=18.6  °C,  Pccalc=54.5  atm;  Tc=9.3  °C, 
Pc=49.7  atm). 

The  mixing  rules  recommended  by  Sanchez  and  Lacombe  (1978)  are  used  to  obtain 
the  mixture  characteristic  parameters.  These  mixing  rules  are  listed  in  Appendix  A.  There 
is  only  one  adjustable  parameter  in  these  mixing  rules,  812,  which  is  directly  fitted  to 
sorption  data.  A  small  value  of  S12  is  used  to  correct  the  deviation  of  the  characteristic 
pressure  of  the  mixture,  Pl2*>  from  the  geometric  mean.  All  of  the  other  parameters  in  the 
mixing  rules  are  pure  component  parameters. 
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Sorption  isotherms  were  calculated  with  a  given  value  of  812  in  the  following 
manner.  At  a  fixed  temperature,  pressure,  and  initial  estimate  of  the  weight  fraction  of  the 
gas  in  the  polymer,  equation  (1)  was  solved  to  obtain  the  volumes  of  the  pure  gas  phase 
and  the  polymer-rich  liquid  phase.  If  equation  (2)  was  not  satisfied,  then  the  weight 
fraction  was  adjusted  and  the  calculational  procedure  was  repeated.  The  best-fit  value  of 
812  was  determined  by  minimizing  the  residual  error  between  calculated  and  experimental 
weight  fractions  for  a  given  sorption  isotherm. 

Results  and  Discussion: 

Sorption  isotherms  were  calculated  for  the  following  systems:  carbon  dioxide- 
PMMA  (Liau  and  McHugh,  1985),  carbon  dioxide-silicone  rubber  (Fleming  and  Koros, 
1986),  ethylene-LDPE  (Cheng  and  Bonner,  1977),  methane-PIB  (Lundberg,  et  al.,  1969), 
methane-LDPE  (Lundberg,  et  al.,  1969),  methane-HDPE  (Lundberg,  et  al.,  1969),  and 
methane-PS  (Lundberg,  et  al.,  1969).  The  pure-component,  characteristic  parameters  used 
in  this  analysis  are  tabulated  in  Table  1  and  the  mixture  parameter,  812,  regressed  from 
experimental  data  for  each  of  the  gas-polymer  systems  is  shown  in  Table  2. 

Figure  1  shows  the  model  calculations  for  the  carbon  dioxide-PMMA  system  at 
41.8,  58.1,  and  68.0°C  (Liau  and  McHugh,  1985;  Wissinger  and  Paulaites,  1986).  The 
812  parameter  obtained  from  the  best  fit  of  this  data  ranges  from  0.019  to  0.035.  For  this 
system  the  fit  of  the  sorption  data  is  very  dependent  on  the  value  of  512.  If  S12  is  set  equal 
to  zero,  a  poor  fit  of  the  data  is  obtained.  Various  values  of  812  will  give  a  reasonable  fit  of 
the  low  pressure  portion  of  the  sorption  isotherms,  but  not  necessarily  of  the  higher 
pressure  portion  of  the  isotherms.  Therefore,  some  data  are  needed  at  high  pressures  to 
define  an  appropriate  value  for  512.  Nevertheless,  the  Sanchez-Lacombe  model  reproduces 
the  trends  in  the  data  and  it  provides  a  means  for  obtaining  reasonable  estimates  of  the 
solubility  of  carbon  dioxide  in  PMMA  at  high  pressures. 
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Although  Koros  (1985)  reports  sorption  data  for  the  carbon  dioxide-PMMA  system 
at  35°C  and  for  pressures  up  to  20  atm,  this  data  is  not  modeled  well  by  the  Sanchez- 
Lacombe  equation  of  state.  However,  the  failure  of  the  Sanchez-Lacombe  model  at  these 
P-T-x  conditions  is  not  unexpected.  Polymers  below  their  glass  transition  temperature 
contain  additional  unrelaxed  volume,  and  PMMA  is  below  its  effective  glass  transition 
temperature  at  the  P-T-x  conditions  investigated  by  Fleming  and  Koros  (1986).  As  shown 
in  Figure  1  (a),  the  carbon  dioxide-PMMA  data  (Liau  and  McHugh,  1985)  is  modeled 
reasonably  well  at  a  slightly  higher  temperature  (41.8°C)  and  higher  pressures  (up  to  300 
atm).  At  the  P-T-x  conditions  employed  by  Liau  and  McHugh  (1985),  the  carbon  dioxide- 
PMMA  system  is  more  than  likely  in  its  rubbery  region.  Although  the  glass  transition 
temperature  of  PMMA  is  105°C,  Chow  (1980)  has  shown  that,  in  general,  the  glass 
transition  temperature  of  a  polymer  can  be  significantly  depressed  if  sufficient  amounts  of  a 
high-pressure  gas  are  sorbed  in  the  polymer.  Based  on  calculations  with  the  expression 
developed  by  Chow  (1980),  it  is  found  that  almost  all  of  the  data  reported  by  Liau  and 
McHugh  (1985)  are  for  PMMA  at  temperatures  above  its  effective  Tg. 

As  mentioned  before,  PMMA  swells  when  in  the  presence  of  high-pressure  carbon 
dioxide.  Figure  2  shows  the  model  calculations  of  this  behavior.  In  these  calculations,  the 
812  values  obtained  from  the  sorption  calculations  are  used.  A  reasonable  fit  of  the 
experimental  data  reported  by  Wissinger  and  Pauliatis  (1986)  is  obtained.  The  model 
predicts  that  if  the  data  of  Wissinger  and  Paulaitis  were  extended  to  higher  pressures  the 
swelling  would  level  off.  This  trend  is  seen  in  the  data  reported  by  Liau  and  McHugh 
(1985),  although  the  absolute  value  of  their  swelling  data  differs  considerably  from  the  data 
of  Wissinger  and  Pauliatus  (1986)  and  from  the  calculated  results. 

Figure  3  shows  the  carbon  dioxide-silicone  rubber  data  reported  by  Fleming  and 
Koros  (1986).  The  812  parameter  determined  for  this  system  is  0.125.  Fleming  and  Koros 
indicate  that  the  silicone  rubber  used  in  their  study  has  a  small  amount  of  cross-linking. 
They  assume  that  since  the  cross-linking  was  minimal,  the  parameters  for  uncross-linked 
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polydimethyl  siloxane  could  be  used  in  the  analysis  of  their  data.  Hence,  the  same 
assumption  is  made  in  this  study.  The  molecular  weight  of  the  silicone  rubber  is  not 
reported  by  Fleming  and  Koros,  so  an  arbitrarily  chosen  molecular  weight  of  50,000  is 
used  in  our  analysis.  As  shown  in  Figure  3,  a  good  fit  of  the  sorption  isotherm  is 
obtained.  The  swelling  of  the  silicone  rubber  in  the  presence  of  carbon  dioxide  is  also 
modeled  using  the  value  of  S12  obtained  from  the  sorption  calculations.  Figure  4  shows 
that  a  good  fit  of  the  swelling  of  the  silicone  rubber  is  obtained  with  the  Sanchez-Lacombe 
equation  of  state.  The  calculated  swelling  for  this  system  exhibits  the  same  trends  that  were 
noted  for  the  carbon  dioxide-PMMA  system.  As  the  partial  pressure  of  carbon  dioxide 
increases,  the  swelling  of  silicone  rubber  increases  monotonically  until  it  eventually  levels 
off  at  pressures  greater  than  80  atm. 

The  fit  of  the  model  to  ethylene-LDPE  isotherms  at  126  °C,  140  °C,  and  155  °C 
(Cheng  and  Bonner,  1977)  is  shown  in  Figure  5.  Like  the  carbon  dioxide-PMMA  system 
presented  earlier,  the  three  isotherms  of  the  ethylene-LDPE  system  are  best  modeled  with 
slightly  different  values  of  812  (612=  0.026, 0.030,  0.041,  respectively).  A  reasonable  fit 
of  all  three  isotherms  is  obtained  as  shown  in  Figure  5.  As  described  earlier,  the  pure 
component  parameters  used  for  ethylene  are  obtained  from  fitting  properties  of  ethylene  in 
the  P-T  space  of  the  experimental  conditions.  If  the  pure  component  parameters  fit  to  the 
vapor-liquid  equilibrium  curve  of  ethylene  are  used  to  calculate  sorption  isotherms,  a 
negative  value  of  812  is  needed  to  fit  the  sorption  isotherms. 

Data  on  the  solubility  of  methane  in  polystyrene  and  polyisobutylene  are  reported 
by  Lundberg,  et  al.  (1969).  Figure  6  shows  a  comparison  of  experimental  and  calculated 
isotherms  for  the  methane-polystyrene  system  at  126  °C,  155°C,  and  188  °C.  The  best  fit 
of  the  data  is  obtained  by  allowing  512  to  vary  as  a  function  of  temperature.  Figure  7 
shows  a  comparison  of  the  methane-polyisobutylene  data  at  155.41°C  with  the  calculated 
results.  Notice  that  even  though  the  partial  pressure  of  methane  is  very  high  for  these  two 
systems,  the  solubility  of  methane  in  the  polymer  is  quite  low.  When  the  sorption 


isotherms  for  these  two  methane-polymer  systems  are  calculated  using  pure  component 
methane  parameters  obtained  from  fitting  methane  properties  at  P-T  conditions  near  those 
of  the  experimental  data,  there  is  not  a  significant  difference  in  the  fit  of  the  data  nor  in  the 
value  of  8j2-  Therefore,  the  methane  characteristic  parameters  given  by  Sanchez  and 

Lacombe  (1976)  are  used  in  the  calculations. 

Figure  8  shows  the  model  calculations  for  the  methane-LDPE  system  (Lundberg,  et 
al.,  1969)  at  155.4  °C  and  the  methane-HDPE  system  (Lundberg,  et  al.,  1969)  at  140.03 
°C.  The  values  of  512  obtained  from  the  best  fit  of  the  data  are  0.0462  and  -0.183, 
respectively.  Notice  that  the  S12  value  for  the  methane-HDPE  system  is  very  negative. 
Even  when  the  pure  component  characteristic  parameters  for  methane  are  fit  to  properties  in 
the  P-T  space  of  the  experimental  data,  the  812  value  remained  negative. 

Conclusions: 

The  solubility  of  high-pressure  gases  in  amorphous  polymers  is  becoming 
increasingly  important  as  membrane  separations  are  being  performed  at  higher  pressures. 
Solubility  data  are  also  important  in  polymer  processes  when  polymers  are  used  as 
controlled-release  devices  or  when  polymers  are  stripped  of  low  molecular  weight 
impurities.  In  this  paper,  we  have  shown  that  the  Sanchez-Lacombe  equation  of  state  with 
only  one  mixture  parameter  can  be  used  to  effectively  predict  the  solubility  of  a  high- 
pressure  gas  in  a  solid  amorphous  polymer  or  a  molten  polymer.  Although  the  mixture 
parameter,  812,  must  be  fit  to  experimental  data,  the  calculations  do  decrease  the  amount  of 
experimental  data  that  is  needed  to  define  a  sorption  isotherm.  Compared  to  using  Flory- 
Huggins  theory  to  predict  the  isotherms,  this  method  avoids  the  need  to  calculate 
hypothetical  states  and  it  rigorously  calculates  activities.  Also,  the  volume  dilation  of  the 
polymer  is  easily  calculated. 
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Appendix  A.  Mixing  rules  used  with  the  Sanchez-Lacombe  equation  of  state. 


As  stated  before,  the  mixing  rules  presented  by  Sanchez  and  Lacombe  are  used  in 
this  analysis.  The  characteristic  pressure  for  the  mixture  is  given  as 


P*=<{)iPi*+<{,2P2*.<j,i<{)2AP* 

(Al) 

where  AP*,  a  mixture  parameter,  is  defined  as 

AP*=P1*+P2*-  2Pj2 

(A2) 

PiVKpj}1”  b-812} 

(A3) 

The  mixing  rules  for  r,  v*,  and  T*  are 

r=x1r1  +X2r2 

(A4) 

v*  =  <|)1v1*+(i)2v2* 

(A5) 

T*  {^+v<|>/T2} 

T  -  1  -<M>2X 

{  4>j  +V(J>2  } 

(A6) 

where 

< 

11 

*|  * 

(A7) 

AP*  v* 

X-  kT 

(A8) 

and  xi  is  the  mole  fraction  of  component  i.  A  detailed  discussion  of  these  mixing  rules  is 
presented  by  Sanchez  and  Lacombe  (1978). 
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Nomenclature: 

k  =  Boltzmann's  constant 

m  =  weight  fraction 

P  =  pressure 

P  =  reduced  pressure 

P12*  =  mixture  parameter 

r  =  segment  length 

5,2  =  mixture  parameter 

T  =  temperature 

T  =  reduced  temperature 

v  =  molar  volume 

v  =  reduced  volume 

X  =  mixture  parameter 

x  =  mole  fraction 

Greek  Letters 

p.  =  chemical  potential 

p  =  density 

p  =  reduced  density 

<|>  =  closed-packed  volume  fraction 

Subscripts 

c  =  critical  properties 

i  =  component  i 

1  =  solvent 

2  =  polymer 

Superscripts 

calc  =  calculated  values 

G  =  gas 

o  =  pure  component 

P  =  polymer 

*  =  characteristic  variable 


Table  1.  Characteristic  equation  of  state  parameters  for  the  light  gases  and  the  polymers 
tabulated  by  Sanchez  and  Lacombe  (1978),  unless  indicated  otherwise. 


Substance 

p*(e/cm3,i 

T*no 

P*  fatml 

CH4 

0.500 

224 

2450 

CO2* 

1.389 

346 

3870 

C02** 

1.452 

290 

5880 

0.511 

336 

1870 

PMMA 

1.269 

696 

4960 

PDMS 

1.104 

476 

2980 

LDFE 

0.887 

673 

3540 

HDPE 

0.904 

649 

4190 

PIB 

0.974 

643 

3490 

PS 

1.105 

735 

3520 

*  Parameters  fit  to  the  liquid  molar  volume,  the  isothermal  compressibility,  and  the 
thermal  expansion  coefficient 

**  Parameters  fit  to  vapor-liquid  equilibrium  line. 


Table  2.  Values  of  the  mixture  parameter  812  used  in  sorption  calculations. 


System 

Temperature 

CC) 

Pressure 

(atml 

Polymer 

Mw 

§12 

CXyPMMA 

41.8 

0-250 

60,600 

0.032 

CCtyPMMA 

58.1 

0-250 

60,600 

0.035 

co2/pmma 

68.0 

0-275 

60,600 

0.019 

co2/pdms*~ 

35.0 

0-60 

50,000 

0.125 

C2H4/LDPE 

126.0 

0-70 

249,000 

0.026 

QH4/LDPE 

140.0 

0-70 

249,000 

0.030 

C^ttyLDPE 

155.0 

0-70 

249,000 

0.041 

ClfyPIB 

155.4 

0-350 

1,300,000 

0.135 

CH4/PS 

125.40 

0-700 

500,000 

0.194 

CltyPS 

155.42 

0-700 

500,000 

0.123 

CH4/PS 

188.40 

0-700 

500,000 

0.099 

CH4/LDPE*** 

155.4 

0-700 

100,000 

0.046 

CltyHDPE*** 

140.03 

0-350 

100,000 

-0.183 

***  The  molecular  weight  for  this  system  the  was  assumed. 
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Figure  1.  Comparison  of  experimental  and  calculated  sorption  isotherms  for  the  CO2- 
PMMA  system.  The  calculated  isotherms  at  41.8°C,  58.1°C,  and  68°C  are 
determined  with  812  equal  to  0.032, 0.035,  and  0.019,  respectively. 
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Figure  3.  Comparison  of  experimental  and  calculated  sorption  isotherms  for  the  C02~ 
silicone  rubber  system.  The  value  of  5*2  used  for  this  isotherm  is  0.125. 
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Figure  5.  Comparison  of  experimental  and  calculated  sorption  isotherms  for  the  C2H4  - 
LDPE  system.  The  calculated  isotherms  at  126°C,  140°C,  and  155°C  are 
determined  with  812  equal  to  0.026, 0.030,  and  0.041,  respectively. 
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Figure  6.  Solubility  of  methane  in  PS  at  125.4°C,  155.42°C,  and  188.40°C  The 

calculated  isotherms  are  determined  with  812  equal  to  0.194, 0.123,  and  0.099, 
respectively. 
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Figure  7.  Comparison  of  experimental  and  calculated  sorption  isotherms  for  the  methane- 
PIB  system.  The  calculated  data  is  determined  with  8j2  equal  to  0.135. 
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Figure  8.  Solubility  of  methane  in  LDPE  and  HDPE.  The  calculated  isotherms  are 
determined  with  8*2  equal  to  0.046  and  -0.183,  respectively. 


